Magnetism and superconductivity at LAO/STO-interfaces: 
the role of Ti 3d interface electrons 
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Ferromagnetism and superconductivity are in most cases adverse. However, recent experiments 
reveal that they coexist at interfaces of LaAlOa and SrTiOa . We analyze the magnetic state within 
density functional theory and provide evidence that magnetism is not an intrinsic property of the 
two-dimensional electron liquid at the interface. We demonstrate that the robust ferromagnetic 
state is induced by the oxygen vacancies in SrTiOa- or in the LaA103-layer. This allows for the 
notion that areas with increased density of oxygen vacancies produce ferromagnetic puddles and 
account for the previous observation of a superparamagnetic behavior in the superconducting state. 

PACS numbers: 74.81.-g,74.78.-w,73.20.-r,73.20.Mf 
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The formation of a metallic state at the interface of 
the bulk insulators LaAlC>3 and SrTiOa 1 has become a 
prototype for the reconstruction of electronic states in 
systems with artificially reduced dimensionality. This 
two-dimensional (2D) electronic system is affected by 
sizable electronic correlations which allow characterizing 
the extended interface electronic states as an electron 
liquid 2 '-. The correlations not only induce a supercon- 
ducting stated but also support magnetism^. A novel yet 
unexplained phenomenon is the coexistence of magnetism 
and superconductivity^ in the 2D electron liquid. 

Recent measurements by Dikin et al. 6 demonstrate a 
hysteretic behavior in the field-dependences of magne- 
toresistance and critical temperature which suggests the 
existence of ferromagnetism in the superconducting sam- 
ples of LaA10 3 (LAO) grown on SrTi0 3 (STO). More- 
over, Dikin et al. take the viewpoint that two distinct 
electronic systems are associated with the antagonistic 
superconducting and ferromagnetic properties: the elec- 
trons that are generated by the polar catastrophe mech- 
anism are suggested to be related with magnetism while 
superconductivity is associated with the charge carriers 
induced by the presence of oxygen vacancies. The con- 
cept of two different types of charge carriers which could 
contribute to the interface electrical transport has been 
discussed also in£~— . 

Very recently, Li et alX probed magnetism through 
torque magnetometry which allows detecting directly the 
magnetic moment of the interface in an external magnetic 
field H. They found a strong superparamagnetic torque 
signal in the superconducting state. With the assump- 
tion that the signal originates from the STO layer next 
to the interface, they obtained the magnetic moment M 
of 0.3 — 0.4/xs per unit cell and a collective magnetic 
moment of the superparamagnetic grains of the order of 
1000 fiB- The observation of a superparamagnetic M(H) 
indicates that ferromagnetic grains form even in the su- 



perconducting state. Magnetic oxygen sites at the A1O2- 
surface (cf. Refi^ 2 .) and the build-up of triplet coupling 
of Ti 3d states through the oxygen bonds (or possibily 
vacancies) in the Ti02 interface plane (cf. Ref^) have 
been proposed^ as scenarios for the formation of a mag- 
netic state. 

The interpretation of the experimental results opens a 
compelling question: can the Ti 3d orbitals that were 
identified in the previous band structure calculations 
(see, e.g., RefsJ^r— ) be responsible both for the metal- 
lic and magnetic states coexisting at the same interface? 
In this paper, we present the results of density functional 
studies which support the existence of a robust ferromag- 
netic state at the LAO/ STO interface induced by oxygen 
vacancies. We demonstrate that both the magnetism and 
conductivity occur involving the Ti 3d electrons but the 
magnetism is due to rather confined electrons around O 
vacancies while the conductivity is a result of the 2D 
electron gas caused by electronic reconstruction. We ar- 
gue that this behavior is a prerequisite of co-existence of 
magnetism and superconductivity which are observed at 
low temperatures. 

To explore whether a ferromagnetic state is induced at 
these interfaces, we consider oxygen vacancies as a mech- 
anism responsible for magnetism. We generate a number 
of supercells which consist of two symmetric LAO/STO 
parts, where each part contains a stack of 4 unit-cell (uc) 
thick LAO layers deposited on a STO-slab of a thickness 
varying between 1.5 uc and 6 uc. The interfacial config- 
uration is considered as Ti0 2 /LaO. The LAO-STO-LAO 
parts are separated by a 13 A thick vacuum sheet. Oxy- 
gen vacancies are assumed to lie in the first interfacial 
Ti02-layers or in one of the A102-layers of the LAO- 
film. A cell with an oxygen vacancy in MO2 (M=TL Al) 
is sketched in Fig. QJb). The vacancy is introduced by 
excluding the oxygen atom 0(a/2, 6/2) in the center of 
the M204-plaquette. The location of the vacancy at the 
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FIG. 1: Schematic view of the SrTiOs/LaA103 heterostruc- 
ture. The supercell contains a 4 unit cell thick LaA103 layer 
deposited on a 2.5 unit cell thick SrTiC>3 slab. The full super- 
cell consists of two symmetric parts of the depicted structure 
and a vacuum layer of 13 A. The structures on the right side 
show (a) a projection of the supercell of STO/LAO on the 
(x,y)-plane of TiC>2, and (b) a M2O4 (M=Ti, Al)-plaquette 
generated for the study of the system with O-vacancies. The 
position of an O-vacancy is identified by a red dashed circle. 



interface is motivated by the experimental evidence of 
O-vacancies present in STO in samples grown at oxygen 
pressures below 10 -5 mba r 23 i 24 . 

Density functional calculations are performed using 
the Generalized Gradient Approximation (GGA) 
in the Perdew-Burke-Ernzerhof pseudopotential 
implementation 25 in the Quantum Espresso (QE) 
package^ 6 -. We use a kinetic energy cutoff of 640 eV and 
the Brillouin zone of the 106- to 166-atom supercells sam- 
pled with 5x5x1 to 9x9x1 fc-point grids. An increase 
of the k-point mesh from (5x5x1) to (7x7x1) leads to a 
negligibly small change of the total energy by 0.005Ry 
and to an increase of the Ti magnetic moments by small 
values of about 0.05 /jb in the presence of O-vacancies. 
The difference between the Ti magnetic moments for 
the two different (2nx2nxl) and (nx2nxl) k-meshes for 
n=4 is about 0.02 /i^, which is a negligibly small value. 
In our calculations we account for a local Coulomb 
repulsion of Ti 3d electrons by employing a GGA+C/ 
approach with Uti — 2 eV— . First, we consider pure 
stoichiometric Ti02/LaO-interfaces as reference for the 
oxygen doped interfaces. The supercells which contain 
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FIG. 2: Projected DOS (in eV" 1 ) for 3d (t 2g ) states of the 
interfacial Ti in the supercell containing 4 unit cell-thick 
LaAlOs layers and a 4-unit cells thick SrTi03 layer. DOS in 
the pure system and in the system with one O-vacancy (25%) 
per supercell area in the interfacial Ti02-layer are shown for 
comparison. The vertical grey line denotes the Fermi level. 



(2x2) planar unit cells have been structurally relaxed 
along the z-axis by a combination of the optimization 
procedures of the full potential WIEN2k-package and 
the pseudopotential QE packag e) 26 ' 27 . The in-plane 
lattice constants have been fixed to their bulk-STO 
cubic values (asTO = &sto = 3.905 A). Similar to 
the previous calculations*^— , we find that a metallic 
state is produced at the LAO/STO interface due to the 
electronic reconstruction. 

Fig. [2] presents the projected Ti 3d densities of states 
(DOS) for both spin directions in the system with su- 
percells containing 4 LAO uc and 4 STO uc along the 
^-direction (the full supercell contains twice the number 
of LAO unit cells). The difference in the spin-projected 
DOS implies a non-zero spin polarization. For the pure 
system without oxygen vacancies, the occupancies of the 
spin-up and spin-down 3d states are almost identical. 
The maximal magnetic moments of the interface Ti are 
about 0.005 fin and the polarization from the more dis- 
tant Ti02 planes is negligible. The calculated magnetic 
moment per (lxl) unit cell of LAO/STO interface is 
0.03 fis which originates mostly from the surface oxygen 
sites. This polarization is too small to support a robust 
ferromagnetic state suggesting that magnetism is not due 
to the pure interface electron gas. 

The situation with O-vacancies is different. An oxygen 
vacancy adds two extra electrons at the interface to pre- 
serve charge neutrality. The two electrons are most likely 
localized in the vicinity of the O-vacancy as found for 
CaO by Elfimov et al. 12 . As shown below, this enhances 
the charge density and increases the exchange splitting 
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FIG. 3: Magnetic moments of Ti atoms in different Ti02 
layers in the (4 uc)SrTiOs/(4 uc)LaA103 structures for the 
structure with one O-vacancy in the interfacial Ti02 layer. 
The black up and red down triangles correspond to the two 
Ti atoms with the planar coordinates (0,0) and (0.5a, 0.5a) in 
a doubled unit cell of SrTiC>3. The TiC>2 layer 4 is the layer 
next to the interface. 



of the spin bands; consequently O- vacancies stabilize the 
ferromagnetic order. 

First, we assume that the oxygen vacancy lies within 
the TiC>2 plane at the interface. In the oxygen-deficient 
system, we find sizable Ti magnetic moments at the in- 
terfacial Ti02-plane (see Fig. The magnetic moment 
of the Ti atoms next to the O-vacancy is ~ 0.33 ps and 
that of the more distant Ti(0,0) atoms at the interfacial 
plane is ~ 0.34 pb- Magnetic moments on Ti atoms away 
from the interface are negligible (Fig. [3|). We also find 
a sizable magnetic moment on the A102 surface plane 
of about —0.18 hb aligned antiparallel to the magnetic 
moment of the interface Ti atoms. Needless to say, the 
concentration of O-vacancies in these model structures is 
far higher than the average density in the experimentally 
investigated heterostructures. Nevertheless, it is evident 
that a triplet coupling is induced between the nearest- 
neighbor Ti sites and that ferromagnetism is enhanced 
in O-vacancy rich regions of the interfacial plane. 

The elimination of the central oxygen in the (2x1) 
configuration results in the formation of stripes of O- 
vacancies along the y-direction characterized by two va- 
cancies near Ti(0.5a, 0) and no vacancies near Ti(0,0) 
atom. To test the stability of such an "inhomoge- 
neous" distribution of O-vacancies, we have also per- 
formed GGA+[7 calculations of a supercell with an or- 
dered 'homogeneous" arrangement of the vacancies cor- 
responding to exactly one vacancy per each Ti atom. 
This can be obtained by elimination of one oxygen in 
the square (y/2 x \/2)-supercell shown in Fig. [TJa). The 
comparison of the calculated total energy with the en- 
ergy for that of the (2 x l)-supercell (Fig. QJb)), both 
containing a 4 u.c. thick-STO layer, gives an energy gain 
of about 0.25 eV per interface u.c. which indicates a ten- 
dency towards an inhomogeneous spatial distribution of 
oxygen vacancies in LAO/STO. 

It is expected that areas with an increased density of 
oxygen vacancies allow for the formation of ferromagnetic 



puddles, as was recently observed by Bert et.alM-. Their 
respective collective magnetic moments are a likely can- 
didate for the source of the superparamagnetic behav- 
ior observed by Li et air . As compared to the ferro- 
magnetically ordered structure, the total energy of the 
(2 x l)-configuration with antiferromagnetically ordered 
interface increases by 0.36 eV/u.c. which implies a high 
stability of the ferromagnetic state. 

The ferromagnetic ordering can be examined in the 
framework of the Stoner model for ferromagnetism. This 
model treats the stabilization of the ferromagnetic state 
as a result of the difference between the reduction of the 
Coulomb interaction for the 3d electrons with parallel 
spins of neighbouring Ti ions and the increase of the ki- 
netic energy caused by the widening of the 3c?-bands for 
the electrons of the same spin^i. The condition for the 
appearance of ferromagnetism is based on the Stoner cri- 
terion Ip(£p) > 1. The interaction / = A/m parame- 
terizes the exchange splitting A, and m is the magnetic 
moment of Ti in the ferromagnetic state (per pb)', p{£f) 
is the paramagnetic density of states at the Fermi level. 
The parameters appearing in the Stoner criterion are de- 
rived from the results of the GGA+ [/-calculations and 
presented in Table HI The Stoner condition is satisfied 
irrespective of the thickness of the STO-layer, which im- 
plies that the ferromagnatic state is favourable in these 
systems although the magnetic moment and exchange 
splitting slightly decrease with increasing STO thickness, 
consistent with Ref. [U The excess 3d charge in the 
vacancy-doped supercells leads to a substantial increase 
of p(ep) up to ~ 1.5 — 1.8 eV -1 as compared to the pure 
systems with p(ef) ~ 0.5 — 0.7 cV _1 . The enhanced 
p{sf) contributes to the strong increase of the Ti mag- 
netic monents in the oxygen-deficient systems. 

O-vacancies strongly influence the electronic structure 
of the Ti 3d states: the excess charge originating from the 
eliminated O atom in the interfacial Ti02 plane leads to 
a redistribution in the occupancy of the five 3c? orbitals. 
The contribution of the 3d e g orbitals to the magnetic 
moment formation is rather insignificant. In contrast, a 
substantial amount of the excess electron charge is trans- 
ferred to the t2 g spin- up orbitals (Fig. [2]). In a (2x l)-unit 
cell, the location of the O-vacancy along the y-direction 
between the two Ti- atoms (see Fig. QJb)) leads to sym- 
metry breaking, splitting the two ti g (3A yz and 3A X2 ) 
orbitals. Due to the random distribution of O-vacancies 



TABLE I: Calculated magnetic moments of interface Ti 
ions nearest to the O-vacancy, the exchange splitting 
(A) and the parameter Ip(ep) of the Stoner criterion for 
(LaA103)4/(SrTi03) n heterostructures with different n. 
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FIG. 4: (a) Magnetic moments of the interface Ti and AIO2- 
layer (per u.c.) in (2STO)/(4LAO)-supercells with a change 
of the position of the O-vacancy along the (OOl)-direction 
between the top AlCh-layer (left) and interface Ti02 (right), 
(b) Variation of the total energy i?tot— -Etot(l) with the change 
of the position of the O-vacancy along the (OOl)-direction. 
Here _E to t(l) is the total energy for the system with an O- 
vacancy located in the top AIO2 surface. 

along the x- and y-directions in LAO/STO samples, the 
electron charge is assumed to occupy both, the 3d yz and 
the 3d xz state — yet the dominant contribution to the 
magnetic moment has to be ascribed to the 3d xy spin- 
up occupancy. 

Oxygen vacancies may appear not only at the interfa- 
cial TiC>2 layer but also on the surface and in the bulk 
of the LAO layer. Consistent with the previous calcula- 
tions — we find that the lowest total energy is achieved 
with the vacancy located in the top A102-surface layer 
(Fig. Hfb)), with an energy gain of about 1.5 cV as com- 
pared to a vacancy in the interface Ti02-layer. The 25 %- 
concentration of surface vacancies produces two electrons 
per (2xl)-unit cell, with one electron transferred to the 
interface and another one hybridized at the AIO2 surface 
close to a central O vacancy, similarly toi£. Our calcu- 
lations show that the placement of an O-vacancy in the 
A102-planes of the LAO-layer still induces a significant 
magnetic moment on the interface Ti atoms (Fig. HJa)). 
This magnetic moment originates from the spin polar- 
ization of the occupied Ti-3d states due to the electron 
charge transferred from the O-vacancy site. Interest- 
ingly, when the O-vacancy lies close to the surface of 
LAO we find a sizable magnetic moment in the AIO2- 



layer with the O-vacancy (Fig. SJa)), similarly to that 
obtained when the O-vacancy is placed in the interfa- 
cial Ti02 layer. This magnetic moment is due to the 
exchange splitting of the not fully occupied local oxygen- 
vacancy-centered state of a mixed Al s-p-character which 
produces quasi-one-dimensional electronic bands and fer- 
romagnetic moments in the A102-layer antiparallel to 
that of the interface Ti atom2£. 

Recent studies of the LAO/STO interfaces distinguish 
two different types of charge carriers in terms of the in- 
terfacial localized Ti 3d electrons and the mobile 3d xy 
electrons of several distant Ti layers^, or they relate 
magnetism to the interfacial electrons produced by the 
electronic reconstruction and associate superconductiv- 
ity to the charge carriers induced by O- vacancies 6 . Our 
findings open a new perspective: both the magnetism 
and the superconductivity are due to the interfacial Ti 
3c? electrons. The magnetism, however, is a result of the 
spin splitting of the populated electronic states induced 
by O-vacancies, while the metallic behavior of the inter- 
face results from the 2D electron liquid caused by the 
electronic reconstruction. The metallic state has been 
related to a superconducting state below 300 mK and 
the predicted scenario suggests that the corresponding 
charge carriers move in regions of small or vanishing O- 
vacancy concentrations. 

Note added during revision: After the submission of 
the manuscript we were informed that Bert et. a/.— 
in fact observed submicrometer puddles of ferromag- 
netism in the presence of the superconducting state at the 
LaA10s/SrTi03 interface. An alternative scenario has 
recently been proposed by K. Michaeli et. alQ who sug- 
gest a homogeneously polarized interface layer which is 
magnetically coupled to a second superconducting layer. 
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